Chiral magnetic effect in lattice QCD with a chiral chemical potential 
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We perform a first lattice QCD simulation including two-flavor dynamical fermion with a chiral 
chemical potential. Because the chiral chemical potential gives rise to no sign problem, we can 
exactly analyze a chirally imbalanced QCD matter by the Monte Carlo simulation. By applying 
an external magnetic field to this system, we obtain a finite induced current along the magnetic 
field, which corresponds to the chiral magnetic effect. The obtained induced current is proportional 
to the magnetic field and to the chiral chemical potential, which is consistent with an analytical 
prediction. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Mh 
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Quantum chromodynamics (QCD) is expected to have 
many characteristic phases at finite baryon density. How- 
ever, it is difficult to extract exact information about 
high-density QCD. One major reason is that lattice QCD 
simulation suffers from the sign problem. At finite baryon 
chemical potential, the fermion determinant becomes 
complex, and its phase fluctuation makes it severe to 
evaluate the ensemble average in the Monte Carlo simula- 
tion. It is instructive to study the exceptional cases which 
can avoid the sign problem in lattice QCD. Through 
these studies, we can learn the qualitative behavior of 
high-density QCD. The well-studied examples are isospin 
chemical potential and two-color QCD. We here consider 
another exceptional case, that is, chiral chemical poten- 
tial [H]. 

The chiral chemical potential \x§ is defined in the 
fermion part of the Euclidean action as 



S F = I d xipdnDp + m + ^57475)^ 



(1) 



This action preserves positive semi-definite property of 
the fermion determinant and thus has no sign prob- 
lem. Physically, the chiral chemical potential generates 
an imbalance between the right-handed and left-handed 
fermions. From the index theorem, this imbalance is 
equivalent to the topological charge of the background 
gauge field. The chiral chemical potential is related to the 
space-time dependent effective theta angle, which repre- 
sents local variation of the topological charge [ll]. The 
constant chiral chemical potential is regarded as a static 
alternative to the topology changing effect. Since lattice 
QCD can simulate only an imaginary-time equilibrium 
system not a real-time evolution, the chiral chemical po- 
tential is a reasonable choice to study this effect in lattice 
QCD. 

One of the most important applications of the chiral 
chemical potential is the chiral magnetic effect. The chi- 
ral magnetic effect is an electromagnetic charged current 
perpendicular to the reaction plane of a heavy-ion colli- 
sion Q . A noncentral collision of two heavy ions produces 
a strong magnetic field, and then this magnetic field fixes 



the spin and momentum directions of the quarks accord- 
ing to their chiralities. As a consequence, a finite net cur- 
rent is induced along the magnetic field if there is a chiral 
imbalance associated with the axial anomaly or, equiv- 
alently, a nontrivial topology of the background gluon 
configuration. The research for the chiral magnetic ef- 
fect is important because it enables us to detect topolog- 
ical structure and local V and CV violation in the strong 
interaction by experiments The chiral chemical po- 
tential has been introduced in several phenomenological 
works for the chiral magnetic effect yj, |j-|6|. The chi- 
ral magnetic effect has also been studied in lattice gauge 
theory without the chiral chemical potential 0-Q ■ 

In this work, we performed a first lattice QCD simu- 
lation with the chiral chemical potential. For the lattice 
gauge action, we used the plaquette gauge action with 
N c — 3. For the lattice fermion action, we used the 
Wilson-Dirac operator as 

D(fJ, 5 )x, y = $x, y - - Jt)U t (x)5 x+ly 

i 

+ (l+H)Ul(x-V)8 x _ l 
- K [(l- 74 e~ )W ^ + ^ 

+(l + 74e -~)C/t (a; _4 )( j x _ a J . (2) 

The chiral chemical potential is introduced as the expo- 
nential matrix factor e ±a A'5T5 _ cosh(a/is) ±75 sinh(a/is), 
which is analogous to baryon chemical potential on the 
lattice [l0(. The lattice action reproduces the original 
action ([1]) in the continuum limit o — > 0. This Dirac 
operator satisfies the relation 75-0(^5)75 = D^/i^) even 
at a 7^ 0, and thus the fermion determinant det D(fi^) 
is always positive real. We generated the gauge config- 
urations with the Nf = 2 dynamical Wilson fermion by 
the Hybrid Monte Carlo algorithm. The inversion of the 
Dirac operator was calculated by the BiCGstab solver 
with the even-odd preconditioning. In most part of the 
following analyses, the lattice gauge coupling is fixed at 
f3 = 2N c /g 2 = 5.32144, and the hopping parameter is 
fixed at k — 0.1665 both for the valence and dynamical 
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FIG. 1: The chiral charge density 715. The lattice sizes are 
xNt = 12 3 x 4, 12 3 x 8, and 12 4 . 



FIG. 2: The expectation value of the Polyakov loop P at 
a/15 = 1.0. The lattice sizes are N! x N t 



8 3 x 4, 12 3 x 4, 



and 16 3 x 4. 



fermions. These values correspond to the lattice spac- 
ing a ~ 0.13 fm (a -1 ~ 1.5 GeV) and the pion mass 
am*- ~ 0.26 (m„ ~ 0.4 GeV) at /x 5 = T = Q. 
First, we calculated the chiral charge density 
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-a 3 (^7475^) = a 3 (ipR>pR 



(3) 



i.e., the difference between the particle number densities 
of the right-handed and left-handed fermions. The nu- 
merical results are shown in Fig. [TJ The chiral charge 
density is exactly zero at fi§ = 0, where the right- 
handed and left-handed fermions are symmetric. The 
chiral charge density becomes nonzero in /zg > and 
increases as 115 increases. Note, however, that the chiral 
charge density saturates around <z//g = 1.4. This is due to 
the saturation of the lattice sites, which is observed also 
in the cases of isospin chemical potential and two-color 
QCD Because the saturation is a lattice artifact, we 
cannot trust the data in a/15 > 1.0. 

Figure Q] includes the data with different values of 
temperature T = l/{aN t ). The low-temperature ones 
(Nt = 8 and 12) are in the confinement phase, and their 
behaviors are almost the same in the present calcula- 
tion. The high-temperature one (N t = 4) is in the de- 
confinement phase. The chiral charge density is qualita- 
tively different between the confinement and deconfinc- 
ment phases. In a/i^ < 0.2, the chiral charge density in 
the deconfinement phase is larger than that in the con- 
finement phase .which is consistent with phenomenolog- 
ical models 0, [l^, 14 1. In a/xg > 0.2, the chiral charge 
density in the confinement phase grows rapidly, and be- 
comes larger than that in the deconfinement phase. This 
rapid growth is considered to be caused by some bound- 
state contribution. 

Next, we focus on the phase structure of a QCD mat- 
ter in ^5 > 0. In two-flavor QCD at 115 = 0, the con- 



finement /deconfinement phase transition is a crossover 
in the non-chiral limit (and a second-order phase transi- 
tion in the chiral limit). There is a possibility that the 
order of the transition is changed by introducing the chi- 
ral chemical potential. For example, phenomenological 
models predict a first-order phase transition line in the 
Hs-T plane 0, Q Ej. 

We calculated the temperature dependence of the 
Polyakov loop, which is a good indicator of confinc- 
ment/deconfinement, by varying the lattice gauge cou- 
pling /3. The numerical result at a/15 = 1.0 is shown in 
Fig. [2] The Polyakov loop increases around j3 ~ 5.27, 
which means a deconfinement transition driven by ther- 
mal effects. However, the Polyakov loop and its sus- 
ceptibility are almost independent of the spatial volume 
V = a 3 N 3 . This scaling behavior indicates that this 
transition is not a true phase transition but a smooth 
crossover. We cannot obtain a signal of a true phase tran- 
sition in a/i5 < 1.0. Therefore, we conclude that, even 
if it exists, there is no change of the order of the phase 
transition in the present accessible range of /15 . This dif- 
ficulty has been already experienced in lattice QCD with 
a finite isospin chemical potential \ii [12| . Although the 
critical endpoint is located away from jij = when the 
quark mass is not quite small. Since the location of the 
critical endpoint depends on the quark mass, the situa- 
tion will change in the near-chiral limit. 

As shown above, we obtained a QCD matter with 
an imbalance between the right-handed and left-handed 
fermion number densities. To analyze the chiral mag- 
netic effect, we applied an external magnetic field to this 
system. On the lattice, the U(l) electromagnetic gauge 
field is introduced as the Abelian phase factor u^[x) on 
the SU(A r c ) link variable U^(x) in the Dirac operator 
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([2]). For generating a constant magnetic field B in the 
X3-direction, the phase factor is taken to be u-i(x) = 
exp(iaqBxi), U\{x\ — aN s ) — exp{—iaqBN s X2), and 
u^{x) = 1 for other components [7[. Only discrete 
values of the magnetic field are allowed as a 2 qB = 
(2n/N^)x (integer) because of the quantization of the 
magnetic flux. Since the magnetic field is not dynamical 
but external, there is no backreaction from the quarks to 
the electromagnetic gauge field. For simplicity, we con- 
sider two fermion flavors with the same charge q. This 
approximation does not change the qualitative behavior 
since the charge difference between the quarks is not es- 
sential for the underlying mechanism of the chiral mag- 
netic effect. 

In this setup, we measured the vector current density 

j M = a 3 (^ 7A1 V) • (4) 

The simulations were done in the deconfinement phase 
(N t = 4), which is relevant for the chiral magnetic ef- 
fect in heavy-ion collisions. The transverse component 
ji and the longitudinal component j'3 of the current den- 
sity are depicted as a function of in Fig. [31 and as a 
function of qB in Fig. |4j The two transverse components 
of the current density are the same, j\ — from the 
rotational symmetry, and they are zero in all of the sim- 
ulations. All components of the current density are zero 
either at B = or at 1x5 = 0. Only when both B and 
/15 are nonzero, a finite current density is generated in 
the longitudinal direction. These results suggest that an 
external magnetic field induces a finite current density 
along the magnetic field only in a chirally imbalanced 
QCD matter. This is exactly what is expected of the 
chiral magnetic effect. 

As seen in Figs. [3] and [4] the induced current density is 
an increasing function of /X5 and qB. Furthermore, it is 
given as a linearly rising function both of /J5 and of qB. 
We can parametrize its functional form as 

j3 = a 3 CN do{ ^ 5 qB . (5) 

The factor A^of is the number of particles with the 
same charge, which is 6 (— N c x Nf) in this simula- 
tion. The overall constant is numerically determined 
as C = 0.013 ± 0.001 by fitting the data. This func- 
tional form has been predicted by an analytical approach 
using the Dirac equation coupled with the background 
magnetic field 1]. The lattice result establishes this 
prediction, except for the overall constant C, which is 
l/(27r 2 ) ~ 0.05 in the analytical approach. This devi- 
ation comes from some QCD corrections. One possible 
candidate is a correction by the renormalization. The 
local vector current fflis not renormalization-group in- 
variant on the lattice [15j . This is very different from the 
vector current in the continuum. Another candidate is 
the dielectric correction, which reduces the induced cur- 
rent [6]. 
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FIG. 3: The transverse current density j\ and the longitudinal 
current density jz as a function of the chiral chemical potential 
/Us. The black dashed line is a linear function ((5]). The lattice 
size is N't X N t = 12 3 x 4. 
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FIG. 4: The transverse current density ji and the longitudinal 
current density j'3 as a function of the magnetic field B. 



The above situation is completely different from the 
standard lattice QCD without the chiral chemical po- 
tential. In the standard lattice QCD, we cannot ob- 
serve the global induced current. Because the current 
itself is zero and only its local fluctuation is nonzero, 
the chiral magnetic effect is studied only through the 
local fluctuation Q- In principle, lattice QCD can repro- 
duce the gauge configuration with a nontrivial topology, 
which gives a finite chiral imbalance via the index the- 
orem. However, the global topological charge (or the 
global chiral charge) per volume is negligibly small, un- 
less one artificially makes the gauge configuration with a 
huge number of topological charge. The magnetic field 
cannot induce the global current in observable amount. 
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magnetic effect is indirectly analyzed through the local 
fluctuation of the current and the local fluctuation is eas- 
ily affected by various contaminations. Although we have 
not addressed the dependence on the quark mass and the 
effect of renormalization, these issues can be studied sys- 
tematically by the standard lattice QCD techniques. It is 
also straightforward to improve details of the numerical 
simulation, such as, lattice action, lattice spacing, lattice 
volume, etc. 

This work was supported in part by the Grant- 
in-Aid for Scientific Research in Japan under Grant 
No. 22340052. The lattice QCD simulations were car- 
ried out on NEC SX-8R in Osaka University. 



FIG. 5: The longitudinal current density j'3 as a function of 
the chiral charge density ng. 



On the other hand, at a finite chiral chemical potential, 
the chiral charge density is finite and independent of the 
volume. Therefore, owing to the introduction of the chi- 
ral chemical potential, we can observe the global current 
induced by the chiral magnetic effect. 

Finally, in Fig. [5j we plot the induced current den- 
sity as a function of the chiral charge density 715 with a 
fixed magnetic field. We can see that the induced current 
density is approximately proportional to n§. In this sim- 
ulation, the magnetic field is very large, qB 3> Under 
the strong magnetic field, the quantum state of a charged 
particle is dominated by the lowest Landau level. While 
the induced current density cannot generally be written 
in a simple analytical function of n§ , the contribution of 
the lowest Landau level can be written as a linear func- 
tion in the chiral limit Q. Although the lattice QCD 
simulation is not in the chiral limit, we can expect that 
the obtained induced current is dominated by the lowest 
Landau level. 

In this study, we have performed the lattice QCD simu- 
lation with the chiral chemical potential, and verified that 
it actually works well. We have succeeded in obtaining 
nonzero current induced by an external magnetic field, 
which is directly related to the chiral magnetic effect. 
This is a great advantage because, otherwise, the chiral 
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